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ABSTRACT 

Polarimetric observations are affected by leakage of unpolarized light into the polarization channels, in a 
way that varies with the angular position of the source relative to the optical axis. The off-axis part of the 
leakage is often corrected by subtracting from each polarization image the product of the unpolarized map and 
a leakage map, but it is seldom realized that heterogeneities in the array shift the loci of the leaked radiation 
in a baseline-dependent fashion. We present here a method to measure and remove the wide-field polarization 
leakage of a heterogeneous array. The process also maps the complex voltage patterns of each antenna, which 
can be used to correct all Stokes parameters for imaging errors due to the primary beams. 
Subject headings: instrumentation: interferometers — instrumentation: polarimeters — techniques: interfero- 
metric — techniques: polarimetric 



1. INTRODUCTION 

The hardware typically used in radio telescopes has the 
great benefit of observing the Stokes Q, U, and V parameters 
simultaneously with Stokes /, but always allows some mixing 
between the polarization channels, as in Figs. 1 and 2. This 
"leakage" is particularly troublesome when it goes from / into 
Q, U, or V, since the polarized signals are usually a small 
fraction of the total intensity, and therefore easily swamped 
by similarly strong leakages from /. 

A radio interferometer uses two or more antennas to mea- 
sure the amplitudes and phases of the electric field imping- 
ing on their receivers. The measurements are stored as "vis- 
ibilities", which are the correlations of the receiver voltages. 
Given some conditions which this article will assume to have 
been met, the visibilities sample the Fourier transform of the 
sky multiplied by the primary beam (directional sensitivity 
function) of the antennas (Clark 1999; Thompson 1999). 

Mathematically, the effect of a pair of antennas A and B on 
the visibilities they observe, VobaAB, is conveniently expressed 
using the Hamaker-Bregman-Sault (Hamaker el al. 1996) for- 
malism, where the four polarizations are combined into a col- 
umn vector The true visibilities are multiplied on the left by 
a set of Jones matrices, each one the outer product of Jones 
matrices for antennas A and B, i.e. 

Dab = Da ® 

represents the on-axis mixing between the nominally orthog- 
onal polarization channels, often called the "D terms". The 
outer product of two matrices M and N, M i8> N, is formed 
by multiplying each entry of M with all of N, and is used, 
along with a complex conjugation of the second factor, to 
bring together the elements from each antenna in a corre- 
lation. Hamaker et al. (1996) explain the algebraic proper- 
ties, including coordinate transformations, of Jones matrices 
and the outer product in more detail. As in Bhatnagar et al. 
(2006), direction dependent effects can also be included, but 
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Fig. 1 . — Conceptual diagram of polarization leakage in an interferometer. 
Each antenna measures two nominally orthogonal polarizations p and q, but 
they are partially mixed before entering the correlator C. 
b: Baseline (separation) between antennas. 
OA: The optical axis (i.e. pointing direction). 
(/, m): Longitudinal and latitudinal offsets peipendicular to OA. 
n: An arbitrary offset in (/, m). 
r^: Voltage pattern of antenna A, factored to exclude polarization 
leakage. 

HPBW: Half Power Beam Width. 

Dg: n independent factor of the polarization leakage of antenna B. 
n dependent factor of the polarization leakage. 
Vobs.AS ■ Observed vector of visibilities in each polarization. 
Although the diagram places and f above the receiver and D below, each 
includes effects from the feed, reflector surface, and receiver support stmts. 
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they must go inside the Fourier integral: 

y°AB = ^AB J ^^ABinWMSIs (n + «,) e'"*"* dn (1) 

where « is a direction on the sky relative to the "phase track- 
ing center", ric- Tic is set electronically, but usually it is cho- 
sen to coincide with the pointing direction of the antennas. S 
is the Stokes matrix, which transforms the sky's Stokes pa- 
rameters. Is - {I, Q, U, V), into the observational polarization 
basis, typically correlations of either circular or linear polar- 
izations. ^AB and ^AB are respectively the wide-field leak- 
age pattern and primary beam for the correlation of antennas 
A and B. They are sometimes multiplied together to form a 
single Jones matrix which is a generalization of the primary 
beam, but the magnitudes of the effects are more easily as- 
sessed if they are kept separate. With the separation, Fab is 
diagonal since it does not mix polarizations in the observa- 
tional basis, and the diagonal elements of H^ab are all one. 

The on-axis portion of the leakage. Dab, is dealt with by 
standard polarimetric calibration techniques, but the leakage 
varies with direction, growing worse toward the edges of the 
primary beam, as in Fig. 2. This paper is concerned with the 
wide-field polarization leakage, 'Vab, and will assume that 
VobsAB has already been corrected by multiplication with D^^. 

If all of the antennas in an array are identical, the effects of 
the primary beam and leakage patterns can be removed in the 
image plane. At the Dominion Radio Astrophysical Observa- 
tory (DRAO) we previously corrected the wide-field contami- 
nation or polarization by multiplying the Stokes / image with 
"leakage maps", and subtracting the results from the mea- 
sured Q and U images, as in Fig. 4. The leakage maps were 
measured by observing the apparent Q/I and U/I of an intrin- 
sically unpolarized source in a grid of offsets from the primary 
beam center (Peracaula 1999). This correction is performed 
completely in the image plane, so we call it the "image-based" 
leakage removal method. It was immediately applicable for 
DRAG'S Synthesis Telescope (ST, Landecker et al. (2000)) 
since its antennas are equatorially mounted and thus its leak- 
age patterns never rotate relative to the sky. The leakage pat- 
terns of an altitude-azimuth mounted telescope such as the 
Very Large Array (VLA) rotate relative to the sky over the 
course of an observation, but the image-based method can still 
be applied if the data are first broken up into a series of snap- 
shots (Cotton 1994). 

Unfortunately, there are differences, known or unknown, 
between the antennas of any real interferometer. The ar- 
ray may be a combination of antennas from originally sep- 
arate telescopes, such as the Combined Array for Research 
in Millimeter- Wave Astronomy (CARMA, Bock (2006)), and 
most very long baseline interferometers. It could also be in a 
transition period where only some antennas have been mod- 
ified, like the partially Enhanced Very Large Array, and/or 
have serious surface errors as at (sub)mm wavelengths. The 
ST is an example of an array where the antennas are simi- 
lar to each other, but with known differences between them. 
The two outermost antennas have 9.14 m diameters with four 
metal struts supporting their receivers, while the other five 
are 8.53 m in diameter with three struts, made of either metal 
or fiberglass. The differences in antenna diameter obviously 
create differences in the half-power beamwidths (HPBWs), 
which at 1420 MHz are 101.8' for the two outer antennas 
and 108.8' for the rest. The variation in the number and 
composition of the struts affects the scattering of incoming 
light, which is an important component of polarization leak- 



age (most of the rest comes from the feeds). 

In polarization images the differences between antennas are 
seen as mismatches between the standard point spread func- 
tion (PSF, or "dirty beam") and the PSF of the leakage. When 
there are phase differences between the leakages of the an- 
tennas, the effective PSF of the leakage is asymmetric (Ek- 
ers 1999) and offset from the peak of the unpolarized emis- 
sion. The effective PSF of the leakage also varies across the 
field, meaning that subtracting a multiplication of the Stokes / 
map with a leakage map cannot fully correct the polarization 
leakage of a heterogeneous array. Additionally, the response 
of each antenna in an array, both in leaked and true radia- 
tion, depends on the scale of the source(s). Resolved features 
have less power at high spatial frequencies, so antennas that 
only participate in long baselines will contribute little leakage 
to them. Unresolved features have no such attenuation with 
baseline length, and elicit an equally weighted mix of leak- 
age from all antennas. Usually leakage maps are measured 
using a bright unresolved object, so in the case of a heteroge- 
neous array their corrections are only accurate for unresolved 
sources.^ Fig. 4 exhibits both of these problems, as can be 
seen in comparison with Fig. 5, the same image corrected with 
the method described in Section 2. The main change for the 
unresolved sources is the presence or lack of surrounding arcs, 
but the supernova remnant also shows a strong difference in 
the on-source residual leakage. 

The work described here aims to improve polarization 
imaging from the DRAO ST. The telescope is engaged in an 
extensive survey of the major constituents of the Interstellar 
Medium, the Canadian Galactic Plane Survey (CGPS, iaylor 
et al. (2003)) which includes imaging in Stokes parameters Q 
and U at 1420 MHz along the plane of the Milky Way. 

We recently measured the real and imaginary parts of the 
leakage patterns for each antenna of the ST (similar to a holo- 
gram measurement, but with finer spacing over a smaller area) 
and have started using them to correct its Q and U observa- 
tions, as seen in Fig. 5. 

2. REMOVING LEAKAGE FROM LINEAR POLARIZATION FOR A 
HETEROGENEOUS ARRAY 

When the polarization leakage (or primary beam) varies 
with both direction and baseline (i.e. antenna pair), there is 
no way to isolate their effects to one of either the image or 
uv planes. Multiplying V'obs.AB on the left by ^^g(n) does not 
work as it does for Dab, because n must be marginalized away 
by integrating with Is . Is is the true intensity distribution of 
the sky, which is unfortunately unknown. A set of Stokes 
/ CLEAN (ilogbom 1974) components makes an acceptable 
substitute, however, both in the replacement of the true sky by 
CLEAN components, and the temporary neglect of Q, U, and 
V. 

Since the correction for a heterogeneous array must be 
added directly to the visibilities, the / model used must match 
the true / visibility function within the sampled part of the uv 
plane. Specifically, it should not be tapered by any sort of 
smoothing in the image plane, and the almost certain discrep- 
ancies between the CLEAN components and true visibility 
function outside the sampled part of the uu plane are imma- 
terial for this purpose. The CLEANed / image should have 
small enough pixels to avoid quantization errors in the com- 
ponent positions, and be CLEANed to at least a moderately 

^ Even then, only if the images are made with the same baseline weighting 
as used for the leakage map. 



Fig. 2. — An example of leakage from Stokes / into Stokes {/ in a CGPS field 
containing the supernova remnant IC443. The thin circle is the 75' radius 
(24% power) cutoff of the usable part of the beam in polarization. The image 
has not been corrected for the sensitivity dropoff of the primary beam, and 
only includes ST data. Note that it has been CLEANed, so the arcs are mostly 
leakage. The grayscale goes from -5 (black) to 5 (white) mjy/beam. 




Fig. 3. — The sources of the / radiation that leaked into Fig. 2. The grayscale 
goes from -12 (white) to 350 (black) mJy/beam. 



faint level. Very faint / emission does not need to be included 
since it will be multiplied by the leakage, typically less than 
a few percent, and it tends to have many more components, 
which would considerably slow down the calculation of the 
correction. Leakage from such emission could be quickly and 
adequately removed by the image-based leakage map method, 
using the CLEAN residual image as the / map. Calculating 
the correction for both Q and t/ of a CGPS field, with a vari- 
able number, on the order of several thousand, of CLEAN 
components, and 1.2 x 10^ visibilities per polarization, takes 
from 15 minutes to overnight on a 2 GHz personal computer. 

Assuming that = (/, 0, 0, 0) in correcting the wide-field 
leakage of Eq. 1 requires some care, since its validity de- 
pends on what Stokes parameters are wanted, and whether 



Fig. 4. — Fig. 2 (same grayscale) after image-based leakage correction. The 
"on-source" correction of unresolved sources is accurate to 1% of /, close to 
the theoretical precision of the measured leakage map, but the arcs around 
strong leakage remain unaffected. Leakage amplitude differences between 
antennas produce rings, and phase differences produce asymmetric arcs. 




Fig. 5. — Fig. 2 (same grayscale) after correcting leakage using measured 
patterns for each antenna. Leakage measurements were made only inside 
the circle, but they have been extrapolated to the edge of the image, which 
works well for clearing up the arcs of sources slightly outside the limit. The 
remaining arcs are primarily due to differences between the primary voltage 
patterns of the antennas. 

the feeds are circularly or linearly polarized. In general each 
measured Stokes parameter is nominally the true Stokes pa- 
rameter, plus first order leakage from two of the other Stokes 
parameters, plus second order leakage from the remaining 
one. This comes from the leakage Jones matrices for each 
antenna having only ones on diagonal, with the leakage terms 
off-diagonal. As a rule of thumb, the true Q and U can be 
thought of as fractions of /, and V as an even smaller frac- 
tion (i.e. second order). With circularly polarized feeds the 
leakage of / into V is second order, and thus possibly of the 
same magnitude as the leakage from linear polarization, but 
the fact that V - (RR - LL)/2 means it is more likely cor- 
rupted by errors in the right and left gains. Linearly polarized 
feeds replace V with one of Q or f/ in a similar situation. If 
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necessary, multiple Stokes parameters can be CLEANed to 
form an estimate of Is, to be iteratively improved using the 
procedure below. 

The visibilities are corrected using the set of Is CLEAN 
components Vc by subtracting 



(2) 



from the visibilities in each polarization at baseline bAsit)- 
VcbjisiOJ is the set of visibilities in each polarization for an- 
tennas A and B at time t for the jth CLEAN component. We 
prefer to use VctAaWj in the form of Stokes parameters in- 
stead of feed correlations since usually only one image (/) 
needs to be CLEANed before applying the correction, ^ab is 
therefore transformed into Stokes form, ^sab- 



4J 



S,AB 



S'H^abS 



(3) 



For the ST, with its circularly polarized feeds, the correction 
is only applied to Stokes Q and U and second order leakages 
are ignored since the leakage from I to Q and U is first order. 
That reduces the used portion of ^s,ab to linear combinations 
of elements of H^a and 4^^, allowing the leakages of / into 
Qor U for a given baseline to be easily calculated on the fly 
from combinations of leakage maps for the individual anten- 
nas instead of storing leakage maps for each combination of 
antennas: 

IpAB — IpA + IpB (4) 

where f is Q or U. Note that the imaginary part would be 
cancelled out if A and B were identical. The Jones matrices of 
individual antennas are in circular coordinates {p - R,q - L), 
so 



Iqa= HJ^,i2 - 4^a,2i, and 

/£/A =-i(4JA,12 + nJ^,2l). 



(5) 
(6) 



These are the leakage patterns that are shown in Figs. 6 to 9. 
Note that the 12 and 21 subscripts refer to the off-diagonal 
elements of the Jones matrix, not baselines between antennas 
1 and 2. 

3. SIMULATED LEAKAGE MAPS 

Ng et al. (2005) calculated theoretical leakage voltage pat- 
terns for the ST's three and four metal strut antennas. Ap- 
plying them to correcting polarization leakage in the CGPS 
(Taylor et al. 2003) confirmed that heterogeneity in the ST 
was having a noticeable effect on the CGPS polarization im- 
ages that was not being corrected by subtracting the Stokes 
/ images multiplied by leakage maps. The correction still 
left significant residuals, however, which was not surprising 
since the simulated patterns were based on an overly simplis- 
tic model of the ST. Some of the three-strut antennas have 
fiberglass supports for their receivers. Treating those as zero 
strut antennas would be incorrect because each receiver box 
has cables running along one of its supporting struts. The 
unknown eff'ective blockage of those cables, along with the 
partial transparency of the fiberglass struts, made measuring 
the actual leakage patterns essential. 

4. ANTENNA PATTERN MEASUREMENTS 

If one antenna. A, in an interferometer points directly at 
a bright isolated source while the others look at it askew, A 
will not have any off-axis leakage or primary beam attenua- 
tion ('+'^(0) = ^a{Q) = 1), and the effective leakage and pri- 
mary beam patterns will be those of the other antennas alone. 



Such offset observations with one antenna on axis are often 
done for hologrammatic measurements of antenna surface er- 
rors, and with two modifications the hologram scheme can be 
adapted to measure the leakage and primary complex voltage 
patterns of each antenna. 

The first modification is to compress the sampling grid of 
offsets. Since there is a Fourier transform relationship be- 
tween the physical features of an antenna and its angular 
power pattern, hologram measurements need to sample a wide 
section of the celestial sphere to resolve small scale errors (i.e. 
a misadjusted panel or smaller) on an antenna. In an antenna 
pattern measurement, however, it is more important to sample 
the main lobe well, so we confined the sampling grid to within 
the first null. In theory^ the antenna patterns should not vary 
any faster with angle than the primary beam. For the ST that 
means its patterns should be fairly smooth on scales smaller 
than approximately a degree, so the measurements were made 
on a grid with 25' spacing out to a maximum distance of 75' 
from the beam center (the extent of beam used for polarization 
mosaics). 

The second modification is only in software, in that the 
antenna patterns come directly from the measured visibili- 
ties, instead of requiring a Fourier transform like surface error 
measurements. The primary voltage pattern of an antenna B 
comes from a observation with an on-axis reference antenna 
A of an unpolarized and unresolved source s: 

y*'^ = (4J^(0) ® H^lin)) (r^(0) ® sis . (7) 

Since the source is effectively /<5(0) the integral of Eq. 1 was 
readily evaluated for Eq. 7. It can be further simplified by not- 
ing that '+'a(0) and Q(0) are identity matrices, and that (un- 
surprisingly, given the physics it represents) the outer prod- 
uct has the redistribution property (Eq. 5 of Hamaker el al. 

(Ma ® Mb) (N^ ® N/j) = (M^Na) ® (MbNb) . 



Eq. 7 becomes: 



(8) 



VTb = (v. ® CVBin) rB{n)v,y} 

r;,„(«) {p,p*s) + r;_22(«) ^B,i2(«) (z'w;;) 
-r* jj(n) 4J* 21 + r* 22(«)</'..?^> 

^B,2i(«)<?..K> + r;,22(«)<?.?'^> 

Vs is {ps,qs), the voltages that s nominally imposes on the 
feeds, s is unpolarized, so {p^q*) - (^sP*) - 0' {PsP*s) - 
{^j^*), reducing Eq. 8 to 



wobs 

'ab - 



(PsP*) + (isIs) 



^*B,u(n) 

-rB.ii(«)n,2i(«) 
r«22(«)^^,i2(«) 



r*22(«) 



The off^-diagonal elements of B's leakage Jones matrix are 
^B,i2(«) = (VfA„/VStJ*, and 



' AB,qpl " AB,qq' 

M^B,21(«) = (VX,/VX/' 



which completely specifies '4'^, since the diagonal elements 
are 1. 

' Both ttie simulations of Ng et al. (2005) and the more intuitive realization 
that objects smaller than the antenna diameter, such as struts, produce features 
broader than the primary beam. 



- 2 




/ {') ' " ' ' ' ' 

Fig. 6. — Real parts of the voltage leakage from / into Q of antennas 1 (top) 
to 7 (bottom) for bands A (left) to D (right). The colorscale goes from -0.05 
(blue) to 0.05 (red). 

Measuring the primary voltage patterns requires knowing 
(PsPI) (- {'i^1*s))- Their diagonal entries (the only nonzero 
ones) can be estimated'* from a regular on-axis observation 
(i.e. (p^(O)p* (0))), so 



' {pAO)p*s(n))^ 
' {qAmql(n))'\ 



, and 



The 1420 MHz feeds of the ST are not offset from the cen- 
tral axes of the antennas, so there should be no difference be- 
tween its \'B,ii{n) and \'B,22in) because of beam squint. We 
therefore collapse its primary voltage patterns from Jones ma- 
trices to a scalar for each antenna: 



goS-axis,BW 



{pA(0)p's(n)) {qAmqlin)) 



{pAiO)pi(0)) {qAmqim) 



/2. 




□□□□ 



To within the noise, since the effects of the primary voltage patterns are 
defined to be whatever is left after on-axis calibration. 
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Fig. 7. — Imaginary parts of the voltage leakage from / into Q of antennas 1 
(top) to 7 (bottom) for bands A (left) to D (right). The colorscale goes from 
-0.05 (blue) to 0.05 (red). 

This approach can even be useful for telescopes with offset 
feeds, such as the VLA, if care is taken to perform all cali- 
bration and self-calibration with / = (pp + qq)/2 instead of 
pp and/or qq individually (conversation with J. Uson, 2006). 
In practice there is some error introduced for wide-field po- 
larimetry by approximating f with a scalar, since although T 
does not mix polarizations in the observational basis, it typi- 
cally does in the Stokes basis. For circularly polarized feeds 
squint mixes / and V for directions away from the pointing 
center. This does not greatly contaminate / since V is almost 
always ~ 0, but is a serious problem for measuring V, espe- 
cially for continuum observations where spectroscopic tech- 
niques cannot help. The ST does have 1-2% leakage from 
/ into V at the half -power level of the primary beam, and al- 
though it could be interpreted as squint the direction of the ap- 
parent squint sweeps through 180° as the frequency goes from 
band A to D. The Robert Byrd Telescope at Green Bank also 
sees a change in the direction of the apparent squint with fre- 
quency (Heiles et al. 2003). Such a variance with frequency is 
inconsistent with the geometrical effect that affects the VLA. 
The ST has only been used to measure V for exceptional cases 
like pulsars and the Sun, that have strong circular polarization. 
Observations that need to measure V off-axis for more weakly 
polarized sources, especially in continuum, will need to apply 
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Fig. 8. — Real parts of the voltage leakage from / into U of antennas 1 (top) 
to 7 (bottom) for bands A (left) to D (right). The colorscale goes from -0.05 
(blue) to 0.05 (red). 

a more extensive treatment. Similarly, when using linearly 
polarized feeds (Saull & Ehle 199t) squint mixes / with Q 
instead of V, making the fnin) = ^iiin) approximation less 
attractive. 

Using g, the primary beam Bsj(n) for a baseline formed by 
correlating antennas s and t is then 

BsAn) = e'off-axis,.5(«)e'off-axis,r(«)- 

Note that the order of s and t matters when antennas s and t 
are not identical. 

Since the patterns are ratios, the requirement above that s 
be unresolved can be loosened to requiring that its size be 
much smaller than the angular scale of variations in the pri- 
mary beam, to avoid smearing the pattern samples. 

With an interferometric array the patterns can be simulta- 
neously measured for all of the antennas except the reference 
antenna (i.e. B is anything but A in the above equations) by 
keeping only the reference antenna pointed at the source while 
the other antennas look at it with the same grid of offsets. The 
patterns of the antenna used as a reference in that set of obser- 
vations can be measured by repeating the observations with a 
different antenna as the reference. 

5. OBSERVED LEAKAGE MAPS 




/ (') 

Fig. 9. — Imaginary parts of the voltage leakage from / into U of antennas 1 
(top) to 7 (bottom) for bands A (left) to D (right). The colorscale goes from 
-0.05 (blue) to 0.05 (red). 

In order to minimize any effects from interference or 
crosstalk the antennas were placed so that the distances be- 
tween them were no smaller than 47 m. Observations were 
made of 3C 147, an unresolved bright source with a flux den- 
sity of 22 Jy at 1420 MHz. 

The beams were sampled on a square grid with 25' spac- 
ing out to a maximum radius of 75' from the beam center. 
The time spent on each spot was varied to achieve approxi- 
mately the same uncertainty for each leakage measurement, 
by making the integration intervals inversely proportional to 
the nominal value of the primary beam: 

_6/2ai-ccos(2-'^^), 



fint(«) oc COS 



HPBW 



(9) 



The on-axis pointing was observed longer because it was ob- 
servationally convenient and it is relatively important since it 
is used to normalize the patterns. 

The entire grid was observed twice, once with antenna 1 as 
the reference antenna, and then again with antenna 7 as the 
reference antenna. That allowed the leakage maps and pri- 
mary voltage patterns of all antennas in the ST to be measured 
without requiring a separate reference antenna. 

The leakage patterns were sampled out to 75' away from 
the beam center, because that is the portion of the beam used 
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by the CGPS. Beyond that limit (the 24% power level of the 
beam) the leakages are expected to be large, and require long 
integration times to measure with the same accuracy. To help 
remove errors that extend within the 75' from objects just out- 
side it, the leakage pattern measurements are extrapolated, us- 
ing a nearest-neighbor method, as far as 120' away from the 
beam center The leakage patterns are also interpolated with 
cubic splines to a grid with 0.20' spacing to match the pix- 
els of the CLEAN component images. An example correction 
with the measured patterns of leakage from / into U is shown 
in Fig. 5. 

6. QUALITY OF LEAKAGE CORRECTION 

Since the form of primary beam used in Equation 9 is not 
necessarily the correct one, the uncertainty in the primary 
voltage pattern for antenna A, g'ofF-axis./i(«), is calculated as: 



1 



'^/,lsamp,A (n)V IcriM 



/a(0) 



(^samps^(O) 
^samps,A('^) 



/off-axis,A 



■(«) 



1/2 



O" /a(0) 

/a(0) ■ 



^off-axis,A(«) gets its name from acting like direction dependent 
factor of A's gain. nsamps,A(«) is the number of samples for 
antenna A in direction n. 

iQAin) is calculated (for an antenna A that comes before the 
reference antenna, B) as 



lQA(n) 



{RA(n)Ll(0)) {LA(n)Rl(0)) 
{LA(n)Ll(0)) {RA(n)Rl(0)) 



{(R + H^AMLmm) ^ {{L + H'AMRmR^m) 



2{L(n)L*(0)) 



2{R(n)R*m 



{R(n)L*m \ l {L(n)R*m \ 
<L(n)L*(0))j„, ^''^ \{R(n)R*ml ^'^ 



Note that the reference antenna is observing on-axis, so it has 
no leakage. The uncertainty in Iqa comes from the noise in 
the receivers: 



S=RR',LL'JiL',LR' 



dl. 



QA 



dS 



-OS 



The source is intrinsically unpolarized, so the crosscorrela- 
tions without leakage, {RL*)„i and {LR*)„i, are zero, and thus 
so are the derivatives of Iqa with respect to RR* and LL* . The 
uncertainty of Iqa reduces to 



{L(n)L*(0)) + {R{n)R*(0)) 

l9'off-axis,A(«)IW 



since antenna A is the off-axis one. cr/^^ has the same form, 
and in our case is identical since ctq = ctu. 

The uncertainties are roughly independent of n because of 
the time weighting, with an average value for antennas 2 to 6 
of 0.0012. The beam centers are an exception, with average 
uncertainties for antennas 2 to 6 of 6 x 10""*. Antennas 1 and 7 
were each used as reference antennas half of the time, so their 
uncertainties are worse by a factor of nearly V2 (ameliorated 
by their slightly larger diameters). 



7. DISCUSSION 

The measured leakage patterns. Figs. 6 to 9, show that al- 
though there is some overall consistency in the patterns, their 
details are unpredictable, both from antenna to antenna and 
from band to band in frequency. Most noticeably, the an- 
tennas with quadrupod receiver supports, 1 and 7, are struc- 
turally nearly identical, but their leakage patterns do not show 
any more similarity to each other than they do to those of 
the tripod antennas. Likely this is because most of the leak- 
age comes not from the struts, but from the feeds. The feeds 
are nominally identical, and their individual flaws are neither 
easily apparent to visual inspection nor tied to the type of an- 
tenna they are mounted on. This suggests that wide-field po- 
larimetry with even nominally homogeneous arrays requires 
measuring the leakage patterns of each antenna, if the needed 
fidelity warrants it. 

Variation of the leakage patterns from band to band is 
prominent in the real parts of the leakage patterns. This rapid 
change with frequency seems surprising at first glance: one 
might expect properties of a waveguide feed to vary quite 
slowly with frequency, and hardly at all across a band that 
is only 2% of the center frequency. The cause appears to be 
the probes used to feed the reflector at 408 MHz; they are 
housed within the 1420 MHz feed (Veidt et al. 1985). Com- 
puted simulations (B.G. Veidt, private communication) indi- 
cate that these probes cause some fine structure in the perfor- 
mance at 1420 MHz. 

The primary voltage patterns. Figs. 10 and 11, reassuringly 
exhibit only the expected dependence on wavelength; namely 
their angular scales are proportional to the observing wave- 
length. Their apparent tight link to antenna structure suggests 
that primary voltage pattern errors are more amenable to cor- 
rection by adjusting the antennas, as is often done using holo- 
grams. Once the primary voltage patterns are known, their 
effect can also be reduced post-observation, even for an in- 
homogeneous array (Bhatnagar et al, 2006). Currently such 
errors are attacked with direction dependent self-calibration 
(modcal, ( yiHis 199Q), also called peeling), which is vulner- 
able to confusing true features on the sky with unwanted ar- 
tifacts. Measuring the antenna patterns with a bright unre- 
solved calibration source instead of through self-calibration 
with a potentially complicated fainter science target removes 
that vulnerability. 

Although we have only tested heterogeneous array leak- 
age correction with equatorially mounted antennas, in prin- 
ciple it would be even easier to adapt it to antennas on 
altitude-azimuth mounts than the image-based leakage map 
method. Since the leakage voltage pattern method already 
deals with visibilities on an individual basis, the only modifi- 
cation needed would be make xj and yj in Equation 2 func- 
tions of time to account for the rotation of the antennas about 
the optical axis relative to the sky as the Earth turns. 

An implicit, but diflicult to avoid, assumption in correct- 
ing for the effect of beam patterns is that the patterns do 
not change with time or observing elevation. The prospect 
of spending observing time on frequent antenna pattern re- 
measurements, possibly for a set of elevations and frequen- 
cies, is unappealing, so there is considerable pressure to en- 
gineer antennas that are stable enough for occasional mea- 
surements to capture most of the effects. The ST antennas 
were not expected to change significantly with time or ob- 
serving direction, but we confirmed their behavior by compar- 
ing recent leakage measurements to the measurements made 
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Fig. 10. — Real parts of the primary voltage patterns of antennas 1 (top) to 7 
(bottom) for bands A (left) to D (right). The grayscale goes from 0.4 (white) 
to 1.0 (black), and the contours go from 0.4 to 0.9 in steps of 0. 1 . 

by Peracaula of the ST's overall leakage amplitude maps at 
21 cm wavelength. There was little change over the inter- 
vening 10 years, despite some surface modifications to a few 
of the antennas. Stability is expected to be a more serious 
problem for larger (as measured in wavelengths) dishes, espe- 
cially if standing waves create a noticeable resonance effect 
in the leakage at certain observing frequencies. Interpolation, 
or theoretical modeling, may be useful for extending the ap- 
plicability of measured maps to additional elevations and/or 
frequencies. Alternatively, if an extremely accurate correc- 
tion is only needed for one bright source within the field of an 
observation, the antenna patterns could be measured at that 
spot immediately before and after the science observation, as 
opposed to mapping the entire main lobe of the antenna pat- 
terns. 



The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
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Fig. 11. — Imaginary parts of the primary voltage patterns of antennas 1 
(top) to 7 (bottom) for bands A (left) to D (right). The colorscale goes from 
-0.1 (blue) to 0.1 (red). The symmetry of antenna 6's patterns suggests that it 
is out of focus. 

agreement by Associated Universities, Incorporated. The Do- 
minion Radio Astrophysical Observatory is operated as a na- 
tional facility by the National Research Council of Canada. 
The Canadian Galactic Plane Survey is a Canadian project 
with international partners. The survey is supported by a 
grant from the Natural Sciences and Engineering Council 
(NSERC). We thank the reviewers for their time and helpful 
comments. 

R. Kothes kindly pointed out IC443 as a good example of 
the new correction method's efficacy, and D. Routledge help- 
fully expanded upon the simulations of Ng et al. We appre- 
ciate the assistance of I.E. Sheehan in facilitating the mea- 
surements and D. Del Rizzo in partially processing the data. 
The processing was also assisted by K. Douglas' program- 
ming and documentation for antenna surface measurements. 
R.R. appreciates J. Uson, W. Cotton, C. Brogan, and D. Balser 
sharing their experience with the VLA and GET. 



REFERENCES 



Bhatnagar, S., Cornwell, T., & Golap, K. 2006, Correction of errors due to 
antenna power patterns during imaging, EVLA Memo 100, National Radio 
Astrophysical Observatory 1 , 7 



Bock, D. C.-J. 2006, in Astronomical Society of the Pacific Conference 
Series, Vol. 356, Revealing the Molecular Universe: One Antenna is Never 
Enough, ed. D. C. Backer, J. M. Moran, & J. L. Turner, 17-+, ADS link 1 



9 



Claik, B. G. 1999, in ASP Conf. Ser. 180: Synthesis Imaging in Radio 
Astronomy II, ed. G. B. Taylor, C. L. Carilli, & R. A. Perley, 1-, ADS 
link 1 

Cotton, W. 1994, Widefield Polarization Correction of VLA Snapshot Images 
at 1 .4 GHz, AIPS Memo 86, National Radio Astrophysical Observatory 1 

Ekers, R. D. 1999, in ASP Conf. Ser. 180: Synthesis Imaging in Radio 
Astronomy II, ed. G. B. Taylor, C. L. Carilli, & R. A. Perley, 321-334, 
ADS link 1 

Hogbom, J. A. 1974, A&AS, 15, 417, ADS link 2 

Hamaker, J. P, Bregman, J. D., & Sault, R. J. 1996, A&AS, 117, 137, ADS 
hnk 1,4 

Heiles, C, Robishaw, T., Troland, T., & Roshi, D. A. 2003, A Preliminary 
Report: Calibrating the GBT at L, C, and X Bands, GBT Commissioning 
Memo 23, NRAO 4 

Landecker, T. L., Dewdney, P. E., Burgess, T. A., Gray, A. D., Higgs, L. A., 
Hoffmann, A. P., Hovey, G. J., Karpa, D. R., Lacey. J. D., Prowse, N., 
Purton, C. R., Roger, R. S., Willis, A. G., Wyslouzil, W., Routledge, D., & 
Vaneldik, J. F. 2000, A&AS, 145. 509, ADS Inik 1 

Ng, T., Landecker, T. L., Cazzolato, F., Routledge, D., Gray, A. D., Reid, 
R. 1., & Veidt, B. G. 2005, Radio Science, 40, 5014 3, 3 



Peracaula, M. 1999, Instrumental polarization coiTections at 1420 MHz 
for the DRAO Synthesis Telescope fields. Tech. rep.. Dominion Radio 
Astrophysical Observatory 1 

Sault, R. J. & Ehle, M. 1996, The ATCA 13-cm Polarimetric Response, ATNF 
Technical Memo 39.3/088, CSIRO 4 

Taylor, A. R., Gibson. S. J., Peracaula, M., Martin, P. G., Landecker, T. L., 
Brant, C. M., Dewdney, P. E., Dougherty, S. M., Gray, A. D., Higgs, L. A., 
Kerton, C. R., Knee, L. B. G., Kothes, R., Purton, C. R., Uyaniker, B., 
Wallace, B. J., Willis, A. G., & Durand, D. 2003, AJ, 125, 1350 1, 3 

Thompson, A. R. 1999, in ASP Conf. Ser 180: Synthesis Imaging in Radio 
Astronomy II. ed. G. B. Taylor, C. L. Carilli. & R. A. Perley, 11-, ADS 
link 1 

Veidt, B. G., Landecker, T. L., Dewdney, P. E., Vaneldik, J. F, & Routledge, 

D. 1985, Radio Science, 20, 1118, ADS link 7 
Wilhs, A. G. 1999, A&AS, 136, 603, ADS link 7 



